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Abstract 

Using a large sample of charm semileptonic decays collected by the FOCUS pho- 
toproduction experiment at Fermilab, we present new measurements of two semilep- 
tonic branching ratios. We obtain values of y(d+^k- J+~n+) = 0-602 ± 0.010 (stat) ± 0.021 (sys) 

and rffit^ffi = °- 54 ± 033 ( stat ) d = °- 048 0ur D+ ~> result 

includes the effects of the s-wave interference discussed in Reference [1] . 
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1 Introduction 



This paper describes new measurements of charm semileptonic decay rates 
into the D + — > K and Df — * final states relative to the hadronic 
decay rates D + — > fT~7r + 7r + and — > 7r + respectively. The D + — > yU + z/ 
width sets the scale of the ^i(O) form factor describing these decays [2] and 
therefore tests Lattice Gauge and quark model calculations. Such tests of the 
overall form factor scale are particularly relevant to the determination of CKM 
matrix element ratios such as | V cd /V cs \ 2 which relates to (for example) T(D + — > 
pfi + u)/r(D + — > K*fi + u). The — > fi + u process has been frequently used 
to estimate the Df branching ratios via assumptions about T(Df — > i^Vi) 
since no high statistics double tag deduced branching fractions for the Df 
presently exist. Frequently the yields are established by counting D+ — > 
7r + decays making a good measurement of T(Df — > fi + v)/T(D+ — > 7r + ) 
important. 

Recently the CLEO Collaboration obtained a new measurement of -t ^t~^K- 
that is somewhat higher than previous measurements. They state in Reference 

[3], that their new value implies an increase in the ratio ^{dZk-i+u) bringing it 
more in line with early quark model estimates and in considerable discrepancy 
with the ISGW2 model [4]. This Letter discusses a more precise determina- 
tion of the y(d+ ^k- 7T + n+) f rom FOCUS and provides a first measurement of 
this ratio that includes the effects of the interfering s-wave described in our 
recent paper [1]. The interference of this s-wave amplitude with the dominant 
D + — > K contribution to D + — > K~n + fi + v significantly distorts the 

angular decay distributions thus affecting the reconstruction efficiency for this 
state. Throughout this paper, unless explicitly stated otherwise, the charge 
conjugate is also implied when a decay mode of a specific charge is stated. 



2 Experimental and analysis details 



The data for this paper were collected in the Wideband photoproduction ex- 
periment FOCUS during the Fermilab 1996-1997 fixed-target run. In FOCUS, 
a forward multi-particle spectrometer is used to measure the interactions of 
high energy photons on a segmented BeO target. The FOCUS detector is a 
large aperture, fixed-target spectrometer with excellent vertexing and particle 
identification. Most of the FOCUS experiment and analysis techniques have 
been described previously [5]. The FOCUS muon system and typical analysis 
cuts are described in References [6] and [1]. 

To isolate the D + —> K^7i + fi + u topology, we required that candidate muon, 
pion, and kaon tracks appeared in a secondary vertex with a confidence level 
exceeding 1%. The muon track, when extrapolated to the shielded arrays of 
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the inner muon system located downstream of both magnets, was required to 
match muon hits with a confidence level exceeding 5%. The kaon was required 
to have a Cerenkov light pattern more consistent with that for a kaon than 
that for a pion by 1 unit of log likelihood, while the pion track was required 
to have a light pattern favoring the pion hypothesis over that for the kaon by 
1 unit [7]. 

To further reduce muon misidentification, a muon candidate was allowed to 
have no or one missing hit in the 6 planes comprising our inner muon system. In 
order to suppress muons from pions and kaons decaying in our spectrometer we 
required that muon candidates have an energy exceeding 10 GeV. Non-charm 
and random combinatoric backgrounds were reduced by requiring both a de- 
tachment between the vertex containing the K~7i + fi + and the primary produc- 
tion vertex of 10 standard deviations and a reconstructed energy between 40 
and 180 GeV. Background from D + — > K~-k + -k + , where a pion is misidentified 
as a muon, was reduced by requiring that the visible mass ttik-k^ < 1-8 GeV/ c 2 . 
In order to suppress background from D* + — > D°tt + — > (K~ fi + u)7i + we re- 
quired M(K-fi + un+) - M(K-/i + u) > 0.18 GeV/c 2 . The momentum of the 
undetected neutrino was estimated from the D + line-of-flight as discussed in 
Reference [1]. 

In addition to these cuts, that we will call "baseline" cuts, we imposed the 
following two cuts on the sample that we use to quote the branching ratio. 
To suppress possible backgrounds from higher multiplicity charm decay, we 
isolated the Kiifi vertex from other tracks in the event (not including tracks 
in the primary vertex) by requiring that the maximum confidence level for 
another track to form a vertex with the candidate be less than 0.1%. To 
suppress background from the re-interaction of particles in the target region 
which can mimic a decay vertex, we required that the charm secondary vertex 
was located at least three measurement standard deviations outside of all solid 
material including our target and target microstrip system. We will call this 
the "out-of-material" cut. 

Apart from the muon cuts and the cut on the M(K~ n + wn + ) — M(K~ 
mass difference, these same cuts were applied to our D + — > K^n + n + sample. 
We required one of the two pions in the D + — > K~tt + tt + final state to be in 
the inner muon system to better match the angular region required for the 
muon in D + — > K^7i + fi + u events. 

We turn next to a discussion of the cuts used in the — > <fi n + u analysis. 
Most of the cuts used for the D + — > K~TT + fi + u and D + — > K~7t + tt + selection 
were used to select our (pjj, + u and 4>ix + samples. Because the Df lifetime is 
shorter than the D + , we reduced our cut on the primary-secondary vertex 
detachment to greater than 5 standard deviations. To further reduce non- 
charm background we required that our primary vertex consisted of at least 
two charged tracks. To further reduce muon contamination to the Df — > fi + u 
state due to decays of pions and kaons in flight, we required that the confidence 
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level that a muon track had a consistent trajectory through the two magnets 
comprising the FOCUS magnetic spectrometer exceeded 5%. 



3 Analysis of the D + — > A n + v final state 

Figure 1 shows the K~n + mass distributions and fits of the signal we obtained 
using two selections of the cuts described above. A very strong A (896) signal 
is present for both samples. To assess the level of non-charm backgrounds, we 
plot the "right-sign" (where the kaon and muon have the opposite charge) 
and "wrong-sign" K~n + mass distributions separately in Figure 1 (a) and (c) 
for the sample with the baseline cuts and the sample with baseline, out-of- 
material, and isolation cuts respectively. We attribute the dramatic decrease 
in the wrong sign component in Figure 1 (c) compared to Figure 1 (a) to a 
reduction in non-charm backgrounds through the use of the out-of-material 
cut. 

The fits to the subtracted K~n + mass distribution shown Figure 1 (b) and 
(d) are used to estimate the yield of D + — > A n + v events. Two subtrac- 
tions are applied. We first subtract the distributions of wrong-sign from right- 
sign events as a means of subtracting non-charm backgrounds that are nearly 
charge symmetric. We then subtract the anticipated right-sign excess predicted 
by our Monte Carlo simulation which simulates all known charm decay back- 
grounds. The Monte Carlo distribution used in this subtraction is scaled by the 
ratio of the fitted yield of D + — > K^n + n + events observed in the data to that 
observed in the Monte Carlo. This second subtraction both simplifies the shape 
of the non-A component of the observed K~n + mass spectrum and corrects 
for the presence of A* 's from charm sources other than D + — > A* /i + z/ such 
as hadronic charm decays where one of the secondaries is misidentified as a 
muon. 

The fits overlayed on Figure 1 (b) and (d) are relativistic Breit-Wigner line 
shapes plus a 1st order polynomial of the form a(rriKn ~ m o) + where m is 
the pole mass of the A*°. The range of the fit is 0.7 < < 1.1 GeV/c 2 . The 
fit parameters were the resonant Breit-Wigner mass, yield, and width and the 
a and (3 background parameters. 

We found the subtracted data better fit a constant width Breit-Wigner with 
a mass and width consistent with the known [8] parameters of the A rather 
than the conventional p-wave form. We performed a variant of the subtracted 
A~7r + fit where the width was allowed to vary with mass according to T(rriKn) = 
T (p/p ) N where p is the kaon momentum in the kaon-pion rest frame at a 
given m Kn , p Q is value of this momentum when the kaon-pion mass equals the 
resonant mass m and the power N was a fit parameter. We found that the 
data was very consistent with the power N = and inconsistent (at the 20 a 
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Fig. 1. We show two D + — » K*\xv signals and fits (a) The K~tt + mass spectra 
subjected to the baseline cuts. The right-sign signal is plotted in a solid line and 
wrong-sign signal is plotted in a dotted line, (b) Fits to the subtracted baseline 
signal (described in the text). Both the wrong-sign, and charm backgrounds are 
subtracted prior to fitting. The fit, described in the text, consists of a Breit-Wigner 
plus a 1st order polynomial. The D + — ► K* /j, + v yield obtained from this fit consists 
of 23 726 events, (c) The -RT - 7r + mass spectra with an additional out-of-material cut 
and secondary isolation cut. The right-sign signal is plotted in a solid line and 
wrong-sign signal is plotted in a dotted line, (d) Fits to the baseline signal with 
additional out-of-material and secondary isolation cut. The D + 
obtained from this fit consists of 11 698 events. 



K*°fj, + u yield 



level) with N = 3 - the value expected for a p-wave Breit Wigner. We plan 
to present, in a subsequent publication, a detailed analysis of the mjf, line 
shape in D + — > K decays including the effects of interference from the 
s-wave amplitude [1] and possible other sources, efficiency variation, and mass 
dependent form factor and barrier corrections. 



Figure 1 (b) shows the presence of a significant negative offset (/?) that is nearly 
absent in 1 (d) which we attribute to a broad structure excess of wrong-sign 
events in our non-charm background. The D + — > K~-k + 'k + yield was deter- 
mined by fitting the K~n + 'K + mass distributions to a Gaussian peak over 
a polynomial background. Figure 2 shows these fits for the baseline sample 
and the sample with the additional out-of-material and isolation cuts. The 
t(d+^k n+u) k raric hing ratio is derived by dividing the fitted K*° yield in 



r(D+^^-7r+7r+) 

D + — > K /i + z/ events by the fitted D + — > K~-k + -k + peak yield in data and 
then comparing to that obtained in our Monte Carlo with a known input 
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Fig. 2. We show two fitted D + — > K _ 7r + 7r + signals and fits (a) The i^ _ 7r + 7r + 
mass spectra subjected to the baseline cuts. The yield obtained from this fit is 
D + yield of 116 748 events, (b) The K^tt + tt + mass spectra with an additional 
out-of-material cut and secondary isolation cut. The yield obtained from this 
fit is D + yield of 65 421 events. The peak near 2.01 GeV/c 2 is from the decay 
D* + — > D°7r + — > (i"C~7r + )7r + where the 7r pion is included in the primary vertex 
thus mimicking the D + — > If ~7r + 7r + decay topology. 



branching ratio. Because our Z) + — > If 7r + yU + z/ simulation model has an s- 



wave amplitude interfering with D + — > K /i + z/, we correct our ratio by a 

factor of 0.945, the fraction of Monte Carlo events due to the D + — > K 

process alone. We compute the numerator of this fraction, by integrating over 

the K~7t + /j + v phase space the model intensity discussed in Reference [1] with 

the s-wave amplitude set to zero. We then divide this numerator by the in- 

— *o 

tensity integral where both the s-wave amplitude and K amplitudes are set 
to their proper values according to Reference [1]. This approach is frequently 
used in charm Dalitz plot analyzes to assess the fractional contributions of 
various quasi-two-body decay channels such as Df — > <J)tt + even though these 
channels should be described by their quantum mechanical amplitudes rather 

than partial decay rates. After this correction, we obtain J+ J+) = 

0.602 ± 0.010 where the quoted error is statistical only based on our samples 
with baseline, out-of-material and isolation cuts. 



Three basic approaches were used to determine the systematic error on y{d+^k- i!+'n+' 
In the first approach, we measured the stability of the branching ratio with 
respect to variations in analysis cuts designed to suppress backgrounds. In 
these studies we varied the secondary isolation cut, the detachment cut, and 
a cut on the number of tracks in our primary vertex. The square root of the 
sample variance of the branching ratio from 8 such cut sets was 0.71 times 
our statistical error. In the second approach, we split our sample according 
to a variety of cuts applied to both the D + — > K~n + fi + b' numerator and 
D + — > K~it + it + denominator and estimated a systematic based on the con- 
sistency of the branching ratio among the split samples. We split our sample 
based on the reconstructed D + momentum, particle versus antiparticle, and 
the sum of the energy of D + secondary tracks that would strike our calorime- 
ter and thus contribute to our hadronic trigger. We chose these variables since 
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they significantly influence our acceptance. For example the D + momentum 
is tied to both our geometrical acceptance, and particle identification efficien- 
cies. The calorimetric energy for a fixed D + momentum is significantly lower 
for the D + — > K~7v + fi + v than the D + — > K~n + n + final state. The maximum 
systematic error for these splits was 1.57 times our statistical error. In the 
third approach we checked the stability of the branching fraction as we varied 
specific parameters in our Monte Carlo model and fitting procedure. These in- 
cluded varying the level of the background Monte Carlo prior to subtracting, 
the power of the momentum dependence of the width in the Breit-Wigner line 
shape as discussed previously, and the value of the three form factor ratios 
that describe the D + — > K decay distribution. We estimate a combined 
systematic from line shape, background level, and form factor systematic that 
is also 1.57 times our statistical error. Combining all three systematic error 
estimates in quadrature we have: 

T{D+ -> K n+v) = Q 6Q2 ± Q ± Q Q21 ^ 



r(L>+ -> K-1T+1T+) 

We believe that our systematic error estimate is conservative. 
4 Analysis of the Df — > final state 



Figure 3 shows the K~K + mass distributions that we obtained using two ver- 
sions of the cuts described above. Figures 3 (a) and (c) compare the K + K~ 
mass spectra for Df — > fi + u candidates and in a charm background Monte 
Carlo normalized to the observed number of Df — > 7r + events where the 
— > fi + u contribution is excluded. The use of the out-of- material and iso- 
lation cuts both significantly reduces the background level and increases the 
agreement between the observed and that predicted background from charm 
sources according to our Monte Carlo. The difference of the fitted yields 
between the data and background Monte Carlo is the yield of that we will 
attribute to —> fj, + u. 

Figures 3 (b) and (d) compare the K + K~ mass spectra for D+ — > ir + can- 
didates in a .D+ signal region and a sideband region based on the K + K~n + 
invariant mass. The difference in fitted yields between the signal and side- 
band K + K~ mass spectra will form our estimate of the — > 7r + yield. 

The branching ratio was derived by dividing the background sub- 

tracted, fitted yield in D+ — > /i + z/ events by the sideband subtracted, fitted 
yield in data and then comparing to the ratio obtained in our Monte Carlo 
with a known input branching ratio. We obtained 1 1^ + ~ > ^, M ,^ = 0.54 ± 0.033 
where the quoted error is statistical only based on our samples with baseline, 
out-of-material and isolation cuts. 
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Fig. 3. We show signals for Df — > <j) l 1+u an d Df — > 7r + (a) The if~.ftT + 
mass spectra in L>+ — » events subjected to the baseline cuts. The 

data is plotted in a solid line and the background Monte Carlo is plotted in 
a dotted line. The difference in the fitted <p yield in the data and that in 
the background Monte Carlo is 2682 events, (b) The K~K + mass spectra in 
Df — ► (p 7r + events subjected to the baseline cuts. The spectra for the Df 
signal region (1.946 GeV/c 2 < tukk-k < 1-996 GeV/c 2 ) is plotted in a solid 
line and the sideband region (1.913 GeV/c 2 < rriRKn < 1-938 GeV/c 2 and 
2.004 GeV/c 2 < m KK7T < 2.029 GeV/c 2 ) is plotted in a dotted line. The differ- 
ence between the fitted </> yields in the signal and sideband region is 4695 events. 

(c) The K~K + mass spectra in Df — > <p p^v events subjected to the baseline, 
out-of-material,and isolation cuts. The background subtracted <p yield is 793 events. 

(d) The K~K + mass spectra in Df — > <j) tt + events subjected to the baseline, 
out-of-material, and isolation cuts. The sideband subtracted yield is 2192 events. 



The systematic error on the t^f~^t^z^ branching fraction was determined 

in a way similar to that used for the D + — > if /x+za The systematic error 
estimate obtained by varying analysis cuts was 1.06 times the statistical er- 
ror. The systematic error estimate by splitting samples was 1.12 times the 
statistical error. The final estimate due to varying fit parameters such as the 
charm background level was 0.94 times the statistical error. Combining all 
three sources we have: 

J = 0.54 ± 0.033 stat ± 0.048 sys 
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5 Summary 



Table 1 _ tQ 

Measurements of the ^m+^R- ^r~T) branching fraction 



Group 


electron 


muon 


This work 










0.602 ± 


0.010 


± 0.021 


CLEO [3] 


0.74 ± 


0.04 


± 


0.05 


0.72 ± 


0.10 ± 


0.06 


CLEO [9] 


0.67 ± 


0.09 


± 


0.07 








E687 [10] 










0.56 ± 


0.04 ± 


0.06 


OMEGA [11] 


0.62 ± 


0.15 


± 


0.09 








ARGUS [12] 


0.55 ± 


0.08 


± 


0.10 








E653 [13] 










0.46 ± 


0.07 ± 


0.08 


E691 [14] 


0.49 ± 


0.04 


± 


0.05 









Table 1 summarizes measurements of the y[d+^k- Sw+) branching fraction for 
electrons and muons. Our measurement is the first one to include the effects 
on the acceptance due to changes in the decay angular distribution brought 
about by the s-wave interference [1] . After correcting the muon numbers by a 
factor of 1.05 to compare to electrons according to the prescription of Refer- 
ence [8] , we find that all values in the table are consistent with their weighted 
average (0.62 ± 0.02) with a confidence level of 19% if systematic errors are 
added in quadrature with statistical errors. Our number is about 1.57 stan- 
dard deviations below the recent CLEO measurement and about 2.1 standard 
deviations above the number obtained by E691 [14]. 

Table 2 

Measurements of the r ^p£*^<p + +*) branching fraction 



Group 


elect 


ron 






muon 


This work 












0.54 ± 


0.033 ± 0.048 


CLE02 [15] 


0.54 


± 


0.05 


± 


0.04 






E687 [16] 












0.58 ± 


0.17 ± 0.07 


ARGUS [12] 


0.57 


± 


0.15 


± 


0.15 






CLEO [17] 


0.49 


± 


0.10 


± 


0.12 







Table 2 summarizes measurements of the ^ * + ~ > , — tt branching fraction. All 
results are remarkably consistent with an average of 0.54 ± 0.04. 
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